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1. Introduction 

The dispersionless KP hierarchy(dKP or Benney moment chain) is 
defined by 

(1) dt„X{p) = {X{p),B4p)}, n = l,2,---., 
where the Lax operator X{p) is 

oo 

(2) A(p)=p + ^^;„+ip-" 

1 

and it can be used to define a set of polynomials: 

5„(p)=[^!M±, t = l,2,3,--- h = x. 

n 

Here denotes a non- negative part of the Laurent series A"(p). 

For example, 

^ p^ ^ p^ 

B2 = — B3 = —+V2P + V3. 

Moreover, the bracket in stands for the natural Poisson bracket on 
the space of functions of the two variables {x,p): 

(3) {f{x,p),g{x,p)} = d^fdpQ - d^gdpf. 
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The compatibility of (0) will imply the zero-curvature equation 

(4) d^rnBn (p) - (p) = {5„ (p) , (p) } . 

If we denote t2 = y and = t, then the equation (jH) for m = 2, = 3 
gives 

V3x = V2y 

VSy = V2t - V2V2X, 

from which the dKP equation is derived (f 2 = v) 

(5) Vyy = {Vt - VVx)^. 

According to the dKP theory [TJ UHl 122 IHZ] , there exists a wave func- 
tion 5'(A, X, t2, ts,- ■ ■) such that p = Sx and satisfies the Hamiltonian- 
Jacobian equation 

dS_ 

dtn 

It can be seen that the compatibility of (jHI) also implies the zero- 
curvature equation (jl)). Now, we expand B„{p) as 



(6) — = i?„(p)|,=5. 



n n ^-^ 

where the coefficients can be calculated by the residue form: 

Gin = -resx=oo{X'Bn{p)dX) = -^reSp=oo(A'+^^^^c?p), 

z + 1 op 

which also shows the symmetry property 

Grin Gni- 

Moreover, from 

df df 

ULn ^*^m 

this implies the integrability of Gin as expressed in terms of the free 
energy JF (dispersionless r function) j37j 



dtidtn ' 



For example, the series inverse to (0) is 

^' ^~ ~~2A2"3A3"4A4 ' 

where J-'in are polynomials of f 2, f 3, ■ ■ ■ , Vn+i and in fact 

, _ ^In _ A" 

iln — TCSp=QQ dp 

n n 
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are the conserved densities for the dKP hierarchy (0). In [HI Ej, it's 
proved that dKP hierarchy (^Q) is equivalent to the dispersionless Hirota 
equation 

(8) Dimx') = -logP^^^-P^^'^ 



where D{X) is the operator Xl^i 



A 

Next, we consider the symmetry constraint Pill] 

n 

(9) = y^^Ci{Si - Sj), 



1=1 



where Si = S{Xi) and Si = S{\i), Aj, Aj are some sets of points, and 
are arbitrary constants. Notice that from ((7j) we know 

On the other hand, by ^ and ©, we also have 

N N 

D{X)T, = D{X)J2c^{S^-S,) = J2c^{DWS^-D{X)Si) 

i=l i=l 

N 

= -} Ci log -, 

~^ p-pi 

where p = p{X),Pi = p{Xi) and pi = p{Xi). 

Hence we obtain the non-algebraic reduction("waterbag" model) |H1 
IT^ of dKP hierarchy 



N 



(10) A = p - log - — ^ = P + Vs+ip ^ 
where 

(11) ^-+1 = 2^ ; — • 

One remarks that in the limit Aj = Aj + Cj, Cj — ^ 0, keeping Qej = 
di to be constant, the Sato function (fTUI) reproduces the Zakharov's 
reduction 



(12) A=p + ^ + ^ + ---+ 



P-Pl P-P2 P-PN 
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which is the algebraic reduction of dKP hierarchy. 
From (|1U|). we have 

1 ^ 

So (ta = y) 

dyX = {X,^Ci{pi - pi)}. 

This will imply 

(13) dypx = d^[^pl + ^Ci{pi-pi)] 

i=l 

1 ^ 

dyPx = d^[-pl + ^Ci{pi-pi)]. 



N 



i=l 



N 



i=l 

For simplicity, in this paper we only consider the case = 1, i.e., 

(14) dypi = 5,[ip? + ci(pi-pi)] 

dyPi = \^p\ + ci (pi - pi )] , 

and the Lax operator ()10|) is truncated to become 

p-pi 



(15) 



A = p — ci log 



P-Pi 



The equation (jl4j) can also be written as the Hamiltonian system 



Ply 




Cl 







- SH3 - 

Spi 


. Ply . 




_ 


_ j_ 

Cl - 




SH3 

- Spi - 



where 5 is the variation derivative and 



^^ = lj dx[ci{pi - Pi) + 3ct(pi - pi, ^ 
A bi-Hamiltonian structure is defined as (for the case of dKP) 



Ply 




" - 




- 5H_ - 




Spi 




Spi 




= Ji 












= h 




. Ply . 




SH3 
- Spi - 




SH 
- Spi - 
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where 

Cl 



■'X 

Cl 



-i 

Cl 

and H is some Hamiltonian such that J2 is also a Hamiltonian operator 
(Jacobi identity), which is compatible with Ji, i.e., Ji + CJ2 is also a 
Hamiltonian one for any complex number c [T^ IT^ [23 12^] • We hope 
to find J2 and the related Hamiltonian H. 

Besides, from flllj] and (fTH|l we know that V2 = v = pi—pi. Then by 
the theory of symmetry constraints of KP 1201 1^ IHS] hierarchy one 
can consider the dispersive corrections for the waterbag model p5|) . 

This paper is organized as follows. In the next section, we con- 
struct the third-order bi-Hamiltonian structure. Section 3 is devoted 
to establishing first-order bi-Hamiotonian structure using WDVV equa- 
tion in topological field theory. In section 4, we discuss the dispersive 
corrections. In the final section, one discusses some problems to be 
investigated. 

2. Third-Order bi- Hamiltonian Structure 

In this section, we investigate the bi-Hamiltonian structure of the 
two component case p4|) . To find the bi-Hamiltonian structure, one 
can introduce the coordinates 

u = pi+pi, v = pi-pi 

to rewrite equation (fT^ as 



(16) 



u 

V 



y 




where 



Hs = dx[ci{pl - Pl) + 3cl{pi - Pif' 

= -j dx[ci ^^ ^ ^ + 3clv'^] = J dxh^ 

and 5 is the variational derivative. One can observe the conserved 
density 

1 Su'^v + ^ „ 2 21 
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has the separable property 



a^/i3 8i%i _ CiV CiV 2^ 1 



where 

^ + 401 4ci 

f V 

Hence we can identify the equation (|TH|l as the generahzed gas dynamic 
Hamiltonian system or separable Hamiltonian system IM] . There- 
fore, according to the separable Hamiltonian theory in j2lE2]; we know 
that the third-order Hamiltonian operator 

J2 = D^U-'D,U-'aD^, 



where 



or 



U.= \ and ' 

Vrr U^, \ 1 



_1 1 / -IJ'{v)v^ 



is compatible with the first order Hamiltonian operator 

/ \ ^ 

J 

So we can write the equation ()16|) as the bi-Hamiltonian structure 

\-tj \it 



Cl 



h5,UUV 

3 

h5,vvu 



where 



= J h^dx = - y reSp=ooi^^dp)dx 

I r 20 

5 / {Clbl - Pi) + lOclipi - Pif + yc2(pi - P^){pI - Pi) 



5 

+ I4ipl-Pl)^}dx 



5 y ^16^ ' ^ 2 ^ 



5 

+ -CiV \dx. 
o 



2 2 2 

■ U V 



WATERBAG MODEL 7 

Next, one will want to find all the conserved density F{u, v) of equa- 
tion (jl6j) . It's not difficult to see that if J F{u,v)dx is a conserved 
quantity if and only if that F{u,v) satisfies the wave equation 

(17) F^., ^™ 



The wave equation (jl7|) leads to two fundamental hierarchies of con- 
served densities jS2lElj 



rjv+i. 



W P. = E (^H-l-2n). '""^'"')''-^ 



n=0 



Here d^^ = dv and ^ acts on all the factors standing to the right 
of it. For example, 

))^ ■ V = dv [/u('y) dv vfi{v)dv]dv. 
Jo Jo Jo Jo 

For reference, we list the first few members of each sequence in the 
generalized gas dynamic system (fTBj): ) = 1 + ^) 

Fq = V, Fi = uv 

F2 = -M^f + -v^ + 2civ'^ 
2 6 

1 1 

F3 = -u^v + u{-v^ + 2civ^) 
6 6 



24 2 ^6 ^^20 9 3 ^ ^ 

and 

Fq = u, Fi = -M^ + -f ^ + 4ci (f In t> — f ) 

F2 = -u^ + u[-v'^ + 4:Ci{v\nv — v)] 
6 2 

-111 12 8 

F3 = —u'^ + -u^[-v^ + 4:Ci{v\nv-v)] + —v^ + -civ^lnv--civ^ 

+ 8clv^\nv ~20clv'^ 
-111 12 8 

120 6 4 ^ ^24 3 9 

+ 8clv^\nv -20clv^] 
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In fact, one can see that 

^^^^ ^ ci(iV - 1)! ^^ ^ ^regp=oo(A^c^p), N>1 

and from |T8|) we also have (for > 1) 

dFN 



du 



(20, ^ . 

Moreover, we notice that the recursion operator 

(21) R = Ja^r' = D^-^DM-' 
is the square of a simpler first-order recursion operator 

R = D,U-\ 

Then we can easily check that , using (fT7|) and ((201), 

dFjL \ ( dFN \ ( dFjL 

R-'^D i = R-'[ e%] =UJ 

dv / \ du / X \ du 

dF^+x 1 = I dF^+i 

du / X V dv 

and similarly for Fat. Hence the dKP hierarchy of ()lfi|l can be obtained 
using the recursion operator, i.e., 

/ t„ \ dv / 

However, the Lax representation of hierarchy generated by F/v is not 
found. 

Also, using the recursion operator ()21|). one can construct a hierar- 
chy of higher order Hamiltonian densities F^, m = 1,2,3,--- , with 
m indicating the order of derivatives on which they depend, and the 
corresponding commuting bi-Hamiltonian system 

(22) I )=^'i 4 j=^^(4 
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where 

Fq = I xvdx 



' ~ 2] ul- ^^iv)v^ 2J ul-{l + '-fW. 

All the flows of (j22p will commute with the generalized gas dynamic 
system For example, for n = 1, we have, after a simple calcula- 

tion, 



R 



1 


fJ.'-^vlvxx+3flulVa:Vxx+IJ.fM'v'^+fl'u^V^-U^Uxx-3lJ.v'^Ua:Ua: 

(n2-AJi;2)3 



where 



dfi 4ci 
^ dv f 2 

Finally, one remarks that there exists a Lagrangian local in the ve- 
locity fields for the equation (fT^ f up to a scaling): 

V^Ut - U^Vt 
L = — 5 5 2v. 

ui - ^ivi 

The local Lagrangian will exist in bi-Hamiltonian structure with a pair 
of first and third order Hamiltonian operators [201 E2| • 

3. Free Energy and Bi-Hamiltonian Structure 

In this section, we investigate the relations between bi-Hamiltonian 
structure and free energy. Then the compatible first-order Hamiltonian 
operators can be constructed. 

Now, we want to find the free energy associated with the dKP hier- 
archy (P) for the Lax operator of the form (fT3|) . Suppose we are given 
two first-order Hamiltonian operators J\ and J2 {d = d^) 

9'\t) g^^t) / I Vf{t) Vf{t) / ^ \ Vfit) Vf{t) 



def 



10 
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They are both Poisson brackets of hydro dynamic type introduced by 
Dubrovin and Novikov fOl • The bi-Hamihonian structure means 
that Ji and J2 have to be compatible, i.e, J = Ji + aJ2 must also 
be a Hamiltonian structure for all value of a. This compatibility 
condition implies that, for any a, the metric is referred as flat pen- 
cil. The geometric setting in which to understand fiat pencil ( or bi- 
Hamiltonian structure of hydrodynamic system) is the Frobenius man- 
ifold [121 113 UH nil • One way to define such manifolds is to construct 
a function F(t^, t^, • • • , f^) such that the associated functions, 

''''^ ~ dpdpdt^' 

satisfy the following conditions. 

• The matrix rjij = Cuj is constant and non-degenerate. This 
together with the inverse matrix t]^^ are used to raise and lower 
indices. On such a manifold one may interpret rjij as a fiat 
metric. 

• The functions c*^ = rj^'^Crjk define an associative commutative 
algebra with a unity element. This defines a Frobenius algebra 
on each tangent space T^M.. 

Equations of associativity give a system of non-linear PDF for F(t) 

i ' C^-^ii J^J-'-v C^-Lrr J^J-rv Oj-'V J- \ ' 



dt^dtPdt^ ' ^t^'^t'y^t" dt^dPdt^ ' ^t^'^tP^t'' ' 

These equations constitute the Witten-Dijkgraaf-Verlinde-Verlinde (or 
WDVV) equations. On such a manifold one may introduce a second 
fiat metric defined by 

(23) g'^ = d'd^W + d^d'F, 
where 

and the contravariant Levi-Civita connection is 

(24) r^^' = d'd^dt,¥, 

This metric, together with the original metric r/*-', define a fiat pencil 
(i.e, 77*-' + ag^^ is fiat for any value of a). Thus, one automatically ob- 
tains a bi-Hamiltonian structure from a Frobenius manifold M.. The 
corresponding Hamiltonian densities are defined recursively by the for- 
mula 

^ = ^ 
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where n > 1, a = 1, 2, ■ ■ ■ , m, and ip^ = rjatt". The integrabihty con- 



ditions for this systems are automatically satisfied when the cf- are 



defined as above. 

For the waterbag hierarchy (fT3j) . it is obvious that 

and those cf^ can be determined by 
where, using (|T9|l . 



on 

^i"^ = = , , ,„ re3p=oo(A"+Mp), n > 0. 
Ci(n + 1)! 



n = 1, c\2 = 4i = 0, 4 = 1 + ^ = /^(^) 



Simple calculations can get 

4ci 

(26) Cj^i = = 0, = = 1) 
By fl2(ij) ■ we can get immediately free energy 

2 o 

(27) + quadratic terms. 

We notice that the free energy ()27|) has no quasi-homogeneity condition 
and, however, the free energy associated with the Benney hierarchy is 
quasi- homogeneous 

After choosing suitable quadratic terms, then from the free energy 
(|77j). using and (j21I), one can construct J2 as follows: 

. 2t2 2t^ 

J2 




2t^ 2t^ 

t'^d^ + dj'^ t^d^ + dj^ 
t'd, + dj' t^d, + dj^ 

Here we remove the non-analytic part of g'^\ i.e., Int^. Also, one can 
verify directly that J2 is a Hamiltonian operator and is compatible with 
Ji. We notice that the constant ci won't appear in J2. 
Now, using the recursion operator 



t^ + dj^d-' e + d^ed-^ 
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one can construct a hierarchy by 

For example, for m = 1 and m = 2, a simple calculation can yield 











\'iJ 


X ^ 







y.2 ) 1*1/ \ 

which are slightly different from the y flow (fT^ and t3(or t) flow of the 
dKP hierarchy respectively: 



t2 




(29) 
where 



1 f, mH^ + it^)H' 



By comparisons between them, one can see that the non-homogeneous 
terms (or higher- order ci terms) of the water-bag hierarchy could be 
removed in the hierarchy ()28|) . In this way, one can say that the hierar- 
chy (PH|) is perturbed, up to some scalings, by the water-bag hierarchy 
with a perturbation parameter Ci. 

Remark: According to the Kodama-Gibbons formulation jTHj, the Rie- 
mann invariants Ai, A2 of (j29j) are given by, 

Ai = A[Ui) = Ci in 



30)\2 = A(u2 = + ci In — =====^— , 

2 ^/{t^y + Acit^ + 

where Ui and U2 are the real roots of ^\p=ui,u2 = 0, i.e., 

"1 = 5 M2 = r , 
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and the characteristic speeds are 

Then the equation (j29|l can be written as 



(31) 



Also, a simple calculation shows that, using (j3n|l . the fiat metric {ds)'^ 
dt^dt^ becomes, in Riemann's invariants, 

{dsf = Tlll{t){d\,f + T^22{t){d\2)\ 

where 

{dpf 1 e 



?7ii(t) = res. 



^22 (t) = res. 



P2 



dp2 |p=«i 

{ dpf _ 1 e_ 

'dT ~ M " 



Since it's known that waterbag reduction (fT3j) is not scaling invariant 
|16j . we can verify that the metric 

^dsf = Hii^dx.Y + '^idx^f 

Ui U2 

(32) = ^(dtl)2 + ^dtldt2^(^+ci)(dt2)2 

is no more flat fH]- Hence from the theory of Darboux-Egrov metric 
[T^ . one believes that there is no first-order bi-Hamiltonian structure 
for (|29|) for However, we know that (j3ip is (semi-)Hamiltonian 

[T^ l^j and it probably will have a compatible non-local Poisson brack- 
ets of hydrodynamic type PH 123 1211 123 I2H] , deserving further inves- 
tigations. 

Finally, one notices that the metric ()32|) can also be obtained using 
the free energy (jTTj) 

W{t\e) = ht'ft' + 2c,{t'flogt' + ht'f + 2c,{t'f. 

Zi o 
4. DISPERSIVE CORRECTIONS 



In this section, one will investigate the dispersive corrections of wa- 
tebag model from the theory of symmetry constraints of KP hierarchy. 
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Let's briefly describe the KP hierarchy [Hj. The Lax operator of KP 
hierarchy is 



n=l 

and the KP hierarchy is determined by the Lax equation ((9„ = Ti 
X) 

dnL = [Bn, L], 

where -B„ = -L" is the differential part of L". For example (T2 = 

y,T3 = T) 

(33) V2Y = ^ + V,x 

(34) V3Y = ^V3XX + V^x + V2V2X 

(35) V2T = ^V2XXX + V3XX + V^x + '2V2V2X- 

Eliminating V3 and V4, we can obtain the KP equation (V2 = V) 

(36) Vt = \vxxx + VVx + S^Vyy, 

which also can be described as the compatibility condition for the eigen- 
f unction 

(37) = {^d'x + Vd + Vs + Vx)(p. 

To get dKP equation one simply takes T„ — > eT„ = t„ in the KP 
equation (jHHjl , with 

9t„ and V{T„) v{tn), 

to obtain the dKP equation when e —>■ 0. Thus the dispersive term 
^Vxxx is removed. Moreover, letting 

(p = exp — 

e 

in we also have the equation (0) for n = 2, 3 

1 

- 2' 

:3 



Sy = ^Sl + V2 
1 ., 

St = -S^ + V2S^ + V3 



when £ — > 0. The compatibility Sty = Syt will yield the dKP equation 
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Since v = pi — pi, from the theory of symmetry constraints of KP 
hierarchy \2(}\ HI] , one can assume the constraint 



(38) V=(ln0i-ln( 



>2 X, 



where 0i and 02 are arbitrary eigenfunctions, i.e., they both satisfy 
equations We remark that if 0i = exp 02 = exp ^ and 

X eX = X, then V{X, Y, T) v{x, y, t) = pi — pi, where 

Pl = Si^, pl=P2 = S2x 

when 6 —>■ 0. 

Then equations (jHTj) become, i = 1,2, 

(39) (piy = ^0ixx + (ln0i - In02)x0i 

(40) (j)iT = ^0ixxx + (ln0i - In02)x0ix 

+ [\/3 + (ln0i -ln02)xx]0i. 
By equations (jH!^ and (jSH)), a simple calculation yields 

(41) ^ [(ln0i)xF-[(ln02)xF 

and by equations (jMj) and one obtains 

(42) y4 = ^{[(ln0i)xf -[(ln02)x]n- 

2 ~2 

One notices that the dispersionless limits of ()41|) and ()42|) are = 

and f4 = ^^^-3^ , respectively. Unfortunately, using (jH^j) and pn|) . after 
a tedious calculation, one has 

[(ln0i)xx]'-[(ln02)xx]' = O, 

which is a contradiction since 0i and 02 are arbitrary eigenfunctions. 
This means that ()4()j) is not a higher-order Lie-Backlund symmetry of 
(jHUj) or the constraint (jHHj) is not admissablej2ni EH]- 
Also, we can consider the following case generalizing ()H8j) 



(43) [l-/(9x)]"Mln0i-ln02 



or 

V-f{dx)V=[\n<Pi-\n<j)2]x, 

where 

f{dx) = o-idx + 0'2dx + ■ ■ ■ + dndx, ca being constants. 
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The dispersionless limit of is also v = pi—pi- Similarly, from (jH!^ 
and IPH) and (jSHI), one can get 

V, = ^[l-/(9xrM[(ln0i)xf -[(ln02)xf} 

n = i[l-/(axrH[(ln0i)xf -[(ln02)xf 
+ 3[f{dx)V']-3V[fidx)V]}. 

Both the dispersionless limits of V3 and V4 are V3 = and ^4 = 

pI-pI 
3 ■ 

Now, using (jH^j) and (pOj) . a lengthy calculation shows that 

(44) 2V[f{dx)Vx] - [f{dx){V')x] = [{ln<P,)xx? - [(ln02)xx]', 

which is also an contradiction. Hence in the general case (j43p . it's also 
not admissible. For example, letting f{dx) = dx, we have (V2 = V) 

(45) V = (l-9J-i[ln0i-ln02]x 

Vs = i[l-9x]-H[(ln0i)xF-[(ln02)xF} 

V, = l[l-dx]-'mn<l>r)xf-[{\n<P2)x? + SVx} 
and ()44|1 becomes 

-2Fi = [(ln</)i)xx]'-[(ln02)xx]' 

or, using (jlSj) . 

n . 1 . , +n ;^ ^ / [(ln02)xx]^ - [(In^Qxx]^ 
(In 01 -ln02)xx = ±(1 - dx)y , 

which puts constraint on 0i and 02- Therefore ()45j) is not admissible. 

5. CONCLUDING REMARKS 

We have investigated the bi-Hamiltonian structure and dispersive 
corrections of the waterbag model for two component. After introduc- 
ing suitable coordinates, one can identify (fTH|) as a separable Hamil- 
tonian system and thus a third-order bi-Hamiltonian structure is ob- 
tained. Also, using the recursion relation of conserved densities , we 
can find the free energy associated with water-bag model in WDVV 
equation of the topological field theory and thus establish the first- 
order bi-Hamiltonian structure. However,the hierarchy constructed by 
recursion operator is not the same as dKP hierarchy. Also, one consid- 
ers the dispersive corrections from the theory of symmetric constraints 
of KP theory. Some calculations show that these dispersive corrections 
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are not admissible. Finally, one remarks that the solutions of the wa- 
terbag model can be found using the hodograph method in [TSl CH] . 

Several questions remain to be overcome. Firstly, from the theory of 
non-local Poisson brackets of hydrodynamic type [221 12H] , one believes 
the bi-Hamiltonian structure of ()29p (or ()16|) Meserves further investi- 
gations, especially that the free energy ^I7\i is not quasi-homogeneous. 
Secondly, as we see in section 4, the integrable dispersive corrections 
of the waterbag model are still unknown. The main difficulty is in the 
quantization of the Lax operator (jlUj) . i.e., p ^ dx- The exact form is 
not clear and needs further investigations [SH] • Ultimately, we hope to 
generalize the results in section 2 to the general case, for example, four- 
component case. But the computation is more involved. One hopes to 
address these questions elsewhere. 
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